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Scale-dependent halo bias due to local primordial non-Gaussianity provides a strong test of single- 
field inflation. While it is universally understood that single-field inflation predicts negligible 
scale-dependent bias compared to current observational uncertainties, there is still disagreement 
on the exact level of scale-dependent bias at a level that could strongly impact inferences made 
from future surveys. In this paper, we clarify this confusion and derive in various ways that there is 
exactly zero scale-dependent bias in single-field inflation. Much of the current confusion follows 
from the fact that single-field inflation does predict a mode coupling of matter perturbations 
at the level of ~ ~ 5 / 3 , which naively would lead to scale-dependent bias. However, we 
show explicitly that this mode coupling cancels out when perturbations are evaluated at a fixed 
physical scale rather than fixed coordinate scale. Furthermore, we show how the absence of 
scale-dependent bias can be derived easily in any gauge. This result can then be incorporated 
into a complete description of the observed galaxy clustering, including the previously studied 
general relativistic terms, which are important at the same level as scale-dependent bias of order 
/nl*"' ~ 1- This description will allow us to draw unbiased conclusions about inflation from future 
galaxy clustering data. 
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1 Introduction 


Primordial non-Gaussianity is a powerful probe of the physics of inflation. In particular, local 
type non-Gaussianity can be used to discriminate between single-field and multi-held models 
of inhation. Single-held models generally predict negligible local non-Gaussianity [1, 2], 7^“' ~ 
(us — 1) <C 1, such that any detection of non-zero local non-Gaussianity would rule out the single- 
held scenario. On the other hand, multi-held models quite generically predict 1/^7'“ | > 1 [3-5], 
setting a natural target for the precision of primordial non-Gaussianity searches at < 1. 

Gurrently, the best constraint on local non-Gaussianity comes from the cosmic microwave 
background (CMB) temperature and polarization bispectra, = 0.8 ± 5.0 (68 % conhdence 
level) [6] , thus showing no evidence for deviations from Gaussianity, while being about an order of 
magnitude away from the order unity target precision. Unfortunately, due to the limited number 
of modes accessible with the CMB, it will not be possible to improve this error bar by more than 
a factor of two with CMB data only [7]. Fortunately, local primordial non-Gaussianity also leads 
to a strong signature in cosmological large-scale structure, in the form of a scale-dependent halo 
bias A5 ng(^) oc k~‘^ [8]. The signal-to-noise of this effect is thus peaked on very large scales 
and therefore is safely distinguishable from the non-linearities of structure formation. While 
systematics in the large-scale galaxy clustering measurement present a challenge, it has been 
shown that large-volume future galaxy surveys can in principle use this scale-dependent bias to 
reach the target precision < 1 [9-14]. 

Since an order unity precision measurement of is thus feasible, testing single-field inflation 
will require an equally precise prediction for scale-dependent bias in that class of models. However, 
currently, there is still disagreement in the literature about whether the scale-dependent bias is 
identically zero, or instead non-zero at the level of an effective non-Gaussianity of order unity. 
The difference in these predictions is small compared to current observational error bars, both 
amounting to the statement that single-field inflation predicts negligible scale-dependent bias, but 
for future surveys it is essential to settle this question. For example, if the theoretical prediction is 
off by order unity, we might end up falsely interpreting a measurement consistent with single-field 
inflation as evidence for multi-field inflation. 

The main focus of this paper is to resolve the confusion about scale-dependent bias in single- 
field inflation and to derive a consistent answer in multiple ways. The main conclusion will be 
that there is exactly zero scale-dependent bias in these models so that any future detection of 
scale-dependent bias would indeed rule out single-field inflation. 

The confusion in the literature derives from the contrast between two types of arguments. 
On the one hand, the single-field consistency conditions are equivalent to stating that, locally, 
a large-scale curvature fluctuation has no physical effect on small-scale physics and is simply 
equivalent to a coordinate transformation. This strongly suggests (correctly) that there can be 
no modulation of the halo number density by the long mode and therefore no scale-dependent 
bias [15-21]. On the other hand, matter perturbations in synchronous gauge, for example, do see 
a mode coupling at the level of an equivalent ~ —5/3 [22], so that the small-scale variance 
is modulated by a large-scale mode, thus suggesting that there is scale-dependent bias [23-28]. 
We will resolve this paradox by noting that the latter is a modulation of the variance of matter 
perturbations on a fixed coordinate scale, but that the local halo number density can only depend 
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on the variance on a given physical scale. We will then explicitly show that when the variance on 
a fixed physical scale is considered, the modulation cancels out, leaving no physical effect, and 
therefore no scale-dependent bias. 

Moreover, we will explain how both the original order unity mode coupling and the later 
cancellation can be easily understood directly in terms of the single-field consistency conditions. 
Specifically, making use to the work of [29-31] , it is easy to implement the initial conditions from 
single-field inflation in any gauge. The long wavelength mode can be introduced by a change of 
coordinates and therefore it is easy to determine its effect on any local quantity knowing only 
how it transforms under a diffeomorphism^. This basic observation has been applied to a number 
of large scale structure observables by a number of authors (see e.g. [34-44]). We apply this logic 
to the small-scale variance of matter perturbations to elegantly derive the mode coupling, but we 
will show that it can even be applied to the halo number density itself to directly obtain the result 
that A6 ng(^) = 0 without referring to the modulation of the variance of matter perturbations. 

In addition to a precise prediction for scale-dependent bias, there are other subtleties in the 
modeling of galaxy clustering, which, if not properly taken into account, could lead to misinter¬ 
pretation of the data. What we observe in a survey are galaxy overdensities based on positions 
derived from observed redshifts and observed angular positions, which is not the same as the 
galaxy overdensity in some arbitrary given gauge. The difference includes redshift-space distor¬ 
tions, but also more subtle relativistic effects that contribute at order unity at the Hubble scale 
and thus are comparable in magnitude to the scale-dependent bias effect from ~ 1. In addi¬ 
tion, one has to be careful that the very definition of bias is gauge independent. These issues have 
all been addressed previously [45-47]. In this paper, we build on these results, in combination 
with the single-field prediction Abj^a{k) = 0, to present a complete prescription for modeling the 
galaxy clustering that can be used to test primordial non-Gaussianity in an unbiased fashion. 

The paper is outlined as follows. In Section 2, we will provide a recipe for the observed galaxy 
overdensity in single-field inflation that includes the proper relativistic definition of galaxy bias, 
the (absence of) single-field scale-dependent bias, and relativistic projection effects. This section 
serves as a summary of our main result and places it in the context of a complete galaxy clustering 
description. In Section 3, we will derive the scale-dependent bias in single-field inflation from 
several perspectives. We will reproduce many results that have been derived in the literate but 
we will show how they can be made compatible, leading to a single, consistent answer. We will 
conclude in Section 4. 

The paper includes two appendices. In Appendix A, we show how to apply the inflationary 
consistency conditions in both Newtonian and synchronous gauge. We explain the slight dif¬ 
ferences that appear from each perspective. In Appendix B, we show explicitly that the mode 
coupling determined by second order perturbation theory is precisely the same as the one deter¬ 
mined by the consistency conditions. In this sense, the results of Section 3 can be derived by a 
brute force calculation. 

Finally, for the remainder of this article, we will drop the “local” superscript on /^l, as we 
will exclusively discuss local type non-Gaussianity. 

^The same results can be derived in the Newtonian limit by [32, 33] using Galilean invariance. 
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2 The Observed Galaxy Bias in Single-Field Inflation 


We wish to unambiguously express, in the context of single-field inflation, the observed galaxy 
density perturbations in terms of underlying (first-order) cosmological perturbations that can 
be computed using Boltzmann codes such as CAMB. This problem can be divided into three 
components. First of all, even if we are given the fluctuations in physical galaxy number density 
in some gauge, 6g, we need to transform this to the observed galaxy density fluctuations, A^, 
which is estimated based on observed redshifts z and angular positions n. Secondly, to obtain 
Sg, we need a way to define galaxy bias that does not suffer from gauge ambiguities. Finally, we 
need to specify the scale-dependence (if any) of the thus defined galaxy bias parameter. 

The main focus of this paper is to compute and clarify the third component, i.e., is there a 
scale-dependent galaxy bias in single-held inhation? We will discuss this in detail in Section 3. 
In the current section, we simply use the main result of Section 3 as part of a complete recipe 
for the observed galaxy clustering in the context of single-held inhation, containing all three of 
the components discussed above. 

First, the relation between observed galaxy overdensity and the density and metric perturba¬ 
tions in a given gauge has been addressed dehnitively in the literature, e.g. [45-52]. For instance, 
in terms of synchronous gauge huctuations, where (in a spatially hat universe) 


ds'^ 


o2(t) 


—+ (1 — 2 ' 0 ) 6^ dx^dx^ + didjx dx^dx^ 


( 2 . 1 ) 


and where velocity perturbations in cold dark matter are set to zero, [46] gives (borrowing here 
some notation from [52]) 


Ag(n, z) = Ag,ioc(n, z) + Ag,K(n, z) + Agj(n, z), 


( 2 . 2 ) 


with the “local” contribution. 


A X 1 

Ag,loc — "9 + 2 


/ 1T dH 2 , f T .2; 

6.-(l+2Q) + —---(l-Q) — 


X + X") 




(2,3) 


and Ag^K and Agj given in [52] as line-of-sight integrals over metric perturbations. We ignore 
the stochastic noise component here. For this expression to be valid, it is crucial that the 
perturbations on the right hand side are evaluated in synchronous gauge (of course, A^ can in 
principle be expressed in terms of perturbations in any gauge). In Eq. (2.3), 6g is the physical 
galaxy overdensity in synchronous gauge, be the evolution bias, be = dln{a^ng)/dlna (with fig 
the mean physical number density), Q is the response to magnihcation dA4, dehned such that 
magnihcation bias causes Ag ^ Ag + QdAi, and D is comoving distance to the galaxy sample. 

Second, we need to relate the galaxy overdensity 5g to the underlying matter (and metric) 
perturbations by a bias prescription. The standard approach is to define the galaxy bias b as 
the proportionality factor between 5g and the matter overdensity <5. However, this makes the 
definition of b gauge-dependent. The natural gauge choice is the synchronous gauge (e.g. [46]), 
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as here slices of constant coordinate time r corresponds to constant proper time so that there is 
no evolution bias. Thus, 

6g = b6 {synchronous gauge) (2-4) 

defines the bias. In an arbitrary gauge, physical density contrast is replaced by a “gauge invariant” 
quantity 

^Invariant _ ^ = h {6 - MHhr) , (2.5) 

^ a In a 

where 5t is the coordinate time shift between synchronous gauge and the gauge under consid¬ 
eration^. Note that the ratio 6g/d is clearly gauge-dependent. In particular, if the true bias 
b as defined above is scale-independent, the ratio 6g/6 in a gauge where constant coordinate 
time slices do not correspond to constant proper time will in general be scale-dependent. It is 
thus important that we have first defined 6 in a gauge-independent manner before we address 
scale-dependent bias. 

The third component that needs to be addressed then is the specification of the galaxy bias b 
defined in Eq. (2.4), and specifically its scale-dependence. To linear order, and for Gaussian initial 
conditions, the bias is scale-independent and can thus be described by a single free parameter, 
b = bi. However, in the presence of primordial non-Gaussianity, there is an additional, scale- 
dependent bias contribution. Specifically, if the non-Gaussianity is of local type, a large-scale 
potential fluctnation will modulate the initial variance of matter fluctuations on small scales, thus 
modulating the number density of halos that these fluctuations will collapse into at a later time. 
Working in synchronons gauge, and going to Fourier space for convenience, 


^ ^ amug 

dg,i = bi5i + - -^ 

a In cT/j 


dipi 




( 2 . 6 ) 


Here, y? is the metric perturbation evaluated in the early nniverse^, related to the matter over¬ 
density S at redshift z by k'^ip = —3HQidm/2T{k)D{z)k‘^ (nsing the Poisson equation), where 
T{k) is the transfer fnnction of matter perturbations, normalized to 1 at low k, D{z) is the linear 
growth function, normalized such that D{z) = l/{l + z) during matter domination, and <7^ is the 
initial variance of small-scale fluctuations smoothed over a comoving length R enclosing the mass 
of the halos of interest (see Eq. 3.14). The quantity Ug is the galaxy number density in a region 
with initial small-scale variance Technically, we have thus assumed here that galaxies occupy 
halos of a single mass. More generally, the dependence on the variance would be replaced 
by a dependence on the variance on a range of scales corresponding to a range of halo masses. 
We have applied an explicit subscript I to stress that we are looking at the effect of large-scale 
fluctuations compared to the fluctuations that collapse into halos. 

Thus, in the presence of mode coupling (non-zero dlnaji/dipi), the bias gets a scale-dependent 
correction, 

b{k) = bi /\b]<jQ{k) , (2.7) 

^i.e. f = T + St, with r the conformal time coordinate in synchronous gauge and f the conformal time in the 
gauge under consideration 

^The description here in terms of ipi is technically only correct if the initial conditions are taken during matter 
domination with the long mode pi outside of the horizon, which is not possible of the mode entered the horizon 
during radiation domination. In Section 3, we will use a more general description in terms of the curvature 
perturbation (i, but here we use pi since this is a very common description in the literature. 


5 









where Abj^a{k) oc k 
parameter, 


If we write the mode coupling in terms of an effective local non-Gaussianity 

( 2 . 8 ) 


dlnaR _ jAb) 
— ^ /nl ? 


difi 


we obtain for the scale-dependent bias^, 


Abj^Q^k) 


O 

^ JNL 


dlnrig 3HQi}rn 

dlnaji 2T{k)D{z)k‘^ 


( 2 . 10 ) 


There has been significant confusion about the level of scale-dependent bias, i.e. of the mode 
coupling characterized by in Eq. (2.10), in single-field inflation. On the one hand, naively, 
is simply equal to the local non-Gaussianity parameter of C. (the comoving curvature per¬ 
turbation on slices of constant density), which in turn is small and given by the single-field 
consistency condition [1, 2], suggesting that the scale-dependent bias is 

if^ = -^K-l), (2.11) 

where is the tilt of the primordial power spectrum. On further inspection, since C and 6 are 
non-linearly related, their non-Gaussianity parameters are not the same. In synchronous gauge, 
it has been found that the non-Gaussianity of matter perturbations is described by an effective 
non-Gaussianity parameter that gets shifted by —5/3 (e.g. [22]). This has led several authors to 
conclude that single-field inflation predicts scale-dependent bias of order unity [23-26] (see also 
[27, 28]), 

if^ = -^-^(ns-l). (2.12) 

On the other hand, however, basic physical arguments suggest that, modulo gradients, the metric 
perturbation ipi cannot have any effect on small-scale physics and its effect is equivalent to a 
coordinate transformation. This strongly suggests that the scale-dependent bias has to be zero 
[15-21]: 

/^f) = 0. (2.13) 

How are these result compatible? The main goal of this paper is to clarify this confusion. In 
the next section, we will calculate A6 ng both using the mode coupling of density fluctuations in 
synchronous gauge and using the argument that corresponds to a simple coordinate transfor¬ 
mation, and will show that both approaches lead to the same result and that the former approach 
can be understood easily in terms of the latter. We aim to present a pedagogical presentation 
in order to settle the confusion as much as possible. The goal of the current section is simply to 
present the end result, which is that in single-field inflation, there is no scale-dependent bias: 


fi^b) _ „ 
Jnl — u 


(2.14) 


^Note that if we assume a universal halo mass function, the response function of the galaxy number density 


simply becomes, 


dlnua 


= (fei - 1) Sc 


din (Jij 

so that Eq. (2.10) reduces to the scale-dependent bias first discovered in [8]. 


(2.9) 
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In general, the observed galaxy clustering is thus fully described by Eqs (2.2) and Eq. (2.3), 
with the galaxy overdensity given by Eq. (2.4), and the bias by Eq. (2.7). In other words, in 
synchronous gauge, the galaxy overdensity entering Eq. (2.3) is 5g(k) = (6i + A5 ng(^)) <J(k) so 
that in single-field inflation. 


A6ng(^) = 0 • 


(2.15) 


Note however that the conversion from physical galaxy overdensity to observed galaxy overdensity 
in Eq. (2.3) does introduce terms that have the same scale-dependence as the scale-dependent 
bias would have. However, if after taking these well understood terms into account, there is 
remaining evidence for a true scale-dependent bias contribution A^nq oc k~‘^, this would rule out 
single-field inflation. 

While we have in this section expressed the observed galaxy overdensity in terms of syn¬ 
chronous gauge variables, our definition of bias was gauge invariant (Eq. (2.5)) so that our result 
can be applied to obtain Ag starting from any gauge. Existing calculations in Newtonian gauge 
can be written in terms of the same bias definition (see e.g. [47]) and therefore no special care is 
needed and we may again use A6 ng = 0. 


3 Single-Field Consistency and Derivations of Scale-Dependent Bias 

In this section, we will provide several derivations of the result A5 ng = 0 for initial conditions 
satisfying the single-field consistency relation. We will start by reviewing the inflationary single- 
field consistency condition itself in Section 3.1. We will then derive the matter density mode 
coupling between long and short wavelengths using the same arguments that led to the inflationary 
consistency condition in Section 3.2 and we confirm the result by explicitly matching the second 
order perturbation theory solution of matter perturbations to the inflationary consistency relation 
(Appendix B). While we consistently find the same non-zero mode coupling in synchronous gauge, 
for example, it does not lead to scale-dependent bias because the modulation of short modes on 
a given physical scale is zero. We will show explicitly that the conversion from short modes on 
a fixed coordinate scale to modes on a fixed physical scale exactly cancels the aforementioned 
mode coupling and our derivation will make it clear that this cancellation is inevitable, as it is 
a direct manifestation of the consistency condition. We will finally see in Section 3.3 how this 
result can be derived more simply by applying the consistency conditions to the galaxies directly. 
All together, we will provide three different perspectives that can be applied in any gauge to 
derive the same physical result. 


3.1 Single-field consistency conditions 


In order to build intuition, let us first review the derivation of the consistency conditions for 
the comoving curvature perturbation on slices of constant energy density [53, 54]. If we work 
on slices of constant density, on super-horizon scales (i.e. neglecting gradients of the curvature 
perturbation), the metric takes the general form. 


ds'^ 





(3.1) 


Dividing the curvature perturbation into a large- and small-scale contribution, C = Cs + Czj the 
single-field consistency conditions are derived by observing that a large-scale perturbation can. 
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to zeroth order in gradients of Q, be generated by a diffeomorphism. Specifically, we can locally 
describe the pertnrbations in a coordinate system (denoted by a tilde) where there is no large-scale 
mode, 

ds2 = g2Chx) ^^ 2) 

where Cs(x), and more generally the local physics, knows nothing of the existence of a large- 
scale mode. We have written only the space-space contribution to the line element (as there are 
only spatial metric perturbations in this gauge), and we do not apply a tilde to r because we 
will not apply any transformation to it. The large-scale mode is then generated by applying a 
diffeomorphism to a new coordinate system (without a tilde), defined by the dilation, 

x = x(l-0), (3.3) 

leaving the time coordinate unchanged. When Q is a constant, this is a pure gauge mode. 
However, to leading order in derivatives we may modify Q —>• Cz(®) and this matches onto a 
physical long wavelength solution up to higher derivative corrections. Since the line element 
itself is a scalar, in the new coordinate system, we then have, to first order in Q, 

ds^ = ds^ = a^ir) dx^ = a^(T) W (1 20) = a^(r) dx^, (3.4) 

where the hnal equality is true to first order^ in Q. 

Comparing the far left-hand side of Eq. (3.4) to the far right hand-side shows first of all that 
the coordinate transformation indeed generates the large-scale curvature perturbation Q, and, 
secondly, that the small-scale curvature perturbation transforms as a scalar, i.e. Cs(x) = Cs(x). 
Hence, in this x coordinate system where the large-scale mode is manifest, 

C(x) = Cs(x(l + Ci)) = C(x)-h0x-VCs(x), (3.5) 

with Cs(x) statistically independent of Q. This is the squeezed limit, single-field consistency re¬ 
lation. Note that, in formulating the mode coupling, it was crucial that the short mode on 
the right-hand side is statistically independent of the long mode. This is equivalent to it be¬ 
ing the small-scale component of a Gaussian held since for Gaussian helds, different scales are 
independent. Thus we could also write Cs ^ Cf- 

We now wish to relate the mode coupling in Eq. (3.5) to a non-Gaussianity parameter /nl- 
The usual ansatz is 

C = C° + ^/»L ((C=f-((C°)t). (3.6) 

where is a Gaussian held. This leads to the mode coupling, 

G = ^ 0 /nl) Cf • (3.7) 

Comparing to the consistency condition, Eq. (3.5), implies that the mode coupling is really 
described by a mode coupling operator, 

/nl = ^ X • V . (3.8) 

D 

®We could have generated the large mode non-perturbatively with the dilation x = x e”""', but since the long 
mode is assumed to be small, we are content to work to first order in 
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However, in a statistical sense, we can replace this operator by an equivalent numerical value. 
The operator x • V corresponds to a rescaling of coordinates so that, statistically, we expect it to 
correspond to the logarithmic derivative of the root-mean-square fluctuation in (g, he. 
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din 



din A: 



(3.9) 


o k^P/-(k) 

(the minus sign follows from going to Fourier space) where A^(A:) = —dimensionless 
power spectrum of C- This is how the consistency relation is commonly phrased. A more rigorous 
derivation of /nl would involve explicitly calculating the effect of the mode coupling /nl on the 
squeezed limit bispectrum, or equivalently, on the modulation of the small-scale variance by the 
long mode, and equating this to the bispectrum predicted by the local ansatz in (3.7). Details 
can be found in many references, including [1, 2]. 

In the following, we will often make use of both forms for the non-Gaussianity, i.e. as an 
operator and as a number, but we emphasize that the latter description is only valid in a statistical 
sense and not at the level of Eq. (3.7). Note also that the statistical equivalence of x • V to 
—dlnVA^(fc)/dln A: is generally valid, but care has to be taken to use the dimensionless power 
spectrum A^(A;) of the field on which the operator acts. For instance, for the matter perturbations, 
we would get x • V —>■ —1/2 (3 -|- n*). This will play an important role in the next subsection. 

The key insight to be learned from this discussion is that, up to derivatives® of Ch can 
understand the effect of the long mode by performing a diffeomorphism. This insight applies 
in any gauge and to any quantity of interest, provided one knows how it transforms. It is this 
feature that is a prediction of single-field inflation which we will call the single-field consistency 
conditions, rather than the specific value of /nl- 

The single-field consistency conditions follow from the fact that during single-field inflation 
there is a single “clock”, see e.g. [2]. This assumes the evolution has reached the inflationary 
attractor solution (although see [56] for a generalization to non-attractor solutions) and ignores 
decaying modes. The consistency conditions remain valid even after inflation, as long as pertur¬ 
bations evolve adiabatically (so that there still is only a single clock) [29, 31]. Adiabaticity in turn 
is maintained as long as gravity obeys the equivalence principle and as long as scales larger than 
the sound horizon are considered for the long mode. Thus, since we assume general relativity, we 
can always apply the consistency conditions while the long mode is super-horizon, and after radi¬ 
ation domination we can even apply them over a large range of scales inside the horizon. We will 
make use of this late-universe application of the consistency conditions in the following section 
to derive the mode coupling of matter perturbations and to derive scale-dependent galaxy/halo 
bias. 

The application to the late universe has recently been highlighted by [31] (based also on [30]), 
where they were applied directly to the matter perturbations in Newtonian gauge. The techniques 
used in large-scale structure follow from an observation of Weinberg [29] (which was generalized 
in [20]) that, in the context of single-field inflation, a long wavelength adiabatic mode can be 
generated by a diffeomorphism that leaves the gauge fixed. We will follow the same logic, although 

®The conformal consistency conditions would allow us to extend this procedure to include one gradient of the 
long mode [20, 55]. We will not include them here in the interest of pedagogy. 


9 





we will specialize to synchronous gauge (see Appendix A for more details on the relation to 
Newtonian gauge). 


3.2 Scale-dependent bias from mode conpling 


We now turn to the calculation of scale-dependent bias in single-held inhation using the mode 
coupling between long and short modes. We again separate perturbations into long and short 
wavelength components, 5 = 6s + stc. We work in synchronous gauge, where the metric is 
given by Eq. (2.1). In this gauge, constant r slices are hypersurfaces of constant proper time. 
Using the fact that galaxies live in dark matter halos, we can then express the local number 
density of galaxies at some time in the late universe in terms of the initial matter perturbations 


as 


Ug = ng{6i,aj{). (3.10) 

In another gauge, different locations at constant coordinate time might be in different stages of 
their evolution so that the galaxy number density would also explicitly depend on the proper time 
ffuctuations. It is thus crucial that we chose a synchronous gauge to start from (see e.g. [46]). 

Eq. (3.10) says, first of all, that Ug has an explicit dependence on the initial 6i. This is because 
the long mode affects the evolution of the small-scale perturbations and therefore the abundance 
of collapsed objects at a later time. Secondly, there is an explicit dependence on the initial 
variance of small-scale fluctuations on a scale corresponding to the halos of interest. In the 
presence of primordial non-Gaussianity, this variance can in principle be modulated by the large- 
scale metric perturbation Q, leading to scale-dependent bias. To quantify the scale-dependent 
bias in the single-field case, we thus need to compute the modulation 


dCi 


(3.11) 


We choose to describe the long mode in terms of the curvature perturbation ^ because it is 
conserved outside of the horizon. It is common to instead write the modulation in terms of 
the initial Newtonian potential during matter domination, ipi = —3/5Q (see also Appendix B). 
However, note that if the long-wavelength mode entered the horizon before matter-radiation 
equality, it is not possible to consider (super-horizon) initial conditions during matter domination. 

The modulation gets two contributions: 


• There is a mode coupling between 6s and the long mode, 

ds = g 0 /nl^ df , (3-12) 

where we have implicitly defined analogously to the non-Gaussianity parameter for 
cf. Eq. (3.7). This gives a modulation of the small-scale variance on a fixed coordinate scale 
R of 

^ ^r ■ (3.13) 

^We are implicitly taking an excursion set formalism approach to halo formation and we are considering galaxies 
living in halos of a fixed mass enclosed in an initial comoving radius R. However, all our conclusions are valid in 
general as long as galaxies live in halos and halos correspond to overdensities in the initial matter distribution. 
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Note that be different than the /nl appearing in Eq. (3.7) because the matter 

density has a nonlinear relation to C- The small-scale variance is 

(^R = j d^ii\W{kR)\‘^ Ps{k) . (3.14) 

Here, W{kR) is the Fourier transform of, e.g., a spherical top-hat function enclosing a 
sphere with radius R and Psik) is the matter power spectrum. 

• The coordinate scale on which the variance is computed should be the radius R corre¬ 
sponding to a fixed physical lenght scale, o(r) R. This means that the coordinate scale R 
is modulated by Q. 

We will now show that the above two effects exactly cancel out in single-field inflation®. 

In synchronous gauge, we can generate the long mode in the same way we did in comoving 
gauge® to derive the inflationary consistency conditions in Section 3.1. We again start in the 
local coordinate frame where the mode Cz is taken out. Here 

p(x,r) = p(r) (^1-h 5s(x,r)^ , (3.15) 

where the perturbation (Is(x, r) is independent of the long mode. Now we again apply the spatial 
coordinate rescaling, Eq. (3.3), that makes the long mode explicit. Since we do not transform r, 
and the matter density is a scalar, we have 

p{t) (1 -h <5(x, r)) = /9(x, r) = p(x, r) = p(r) (l + 5(x, r)^ , (3.16) 

so that 6s transforms as a scalar (as did Cs) under the dilation, 

(5s(x,r) = 5^(x,r). (3.17) 

Thus, in the coordinate system with the long mode explicit, the mode coupling takes the same 
form as for 

^^(x, r) = 5s(x(l -h Cl)) = ds{x) + Cz X • V55 (x) , (3.18) 

or, 

Jp) =2Czx-(3.19) 

Here we have expanded the matter perturbation <5 = -|- ^6^'^'^ -|- .... An alternative way to 

derive this mode coupling, see e.g. [23-26], is to compute the second order perturbation theory 
solution for matter perturbations in the post-inflationary epoch and to match this solution to the 
consistency condition for C given by single-field inflation. We follow this approach for the special 
case of a matter dominated Universe in Appendix B and show that it yields the exact same result 
found in a more straightforward way above. 

®In many cases, the above procedure is performed in Lagrangian coordinates. In that case, it has been shown 
that the consistency conditions acts trivially [41]. 

®The technical subtleties that have been discussed in comoving gauge and Newtonian gauge for obtaining 
physical solutions are not important in synchronous gauge. See Appendix A for details. 
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Comparing this to the definition of the effective non-Gaussianity for the matter perturbation, 
Eq. (3.12), 

^ = “E Ci /nl , (3.20) 

5 

we find = | x-V. Statistically, the operator x-V acting on gives —= —5 i^ + rig) 
(since the matter power spectrum oc oc Therefore, we get 

/nl “ “ x 2 

This is the squeezed limit mode coupling result found in the literature, e.g. [22]. 

This leads to the following modulation of the small-scale variance on a given coordinate scale 

R, 

{'^r)ipi ~ = (1 “ 40 — {us — 1 ) Cl) a% . ( 3 . 22 ) 

However, the halo density should depend on the small-scale perturbation statistics on a fixed 
physical scale. The long-wavelength part of the spatial component of the metric is given by. 


ds'^ 



dx^ -|- didjXi dx'^dx^ 


«^('^) [(1 + 2 0 ) d -^^] , 


(3.23) 


where in the second line we have used Oi = ~0 under the diffeomorphism defined in Eq. 3.3, and 
we have again neglected gradients of the long mode. Thus a physical distance is given in terms 
of coordinate distance in synchronous gauge by 


a(T) R = a{T) (1 -h 0) R- 


(3.24) 


This of course simply reflects the coordinate rescaling Eq. (3.3) we applied to generate the long 
mode in the first place. The small-scale variance on a fixed physical scale a(r) R is now 

2 _ 


Here, as in the above, we treat the power spectrum as a pure power law so that d\nP^{k)/dink = 
Us , but is is trivial to extend this to power spectrum with running of the spectral index and higher 
order corrections. 

The variance on the right-hand-side above, is evaluated on a fixed coordinate scale, but 
still gives the small-scale variance at that scale in the synchronous gauge coordinate system. 
Eq. (3.22) gives the modulation of that variance by Cl- Putting it all together, we find for 


d^k\W{k R{1 + 5/3ipi))\^ Ps{k) 

d^klWikm^PsCkC^ + Ci)) 


a 


R 


1-1-30 + 0 


din Ps{k) 
dink 


[1 + 40+ (us — 1) 0] • 


(3.25) 
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the variance on a fixed, physical scale, and explicitly writing the dependence of the small-scale 
variance on Q, 

= [1 + ‘‘Ci + (n. - 1) C.1 (4)^_ 

= [1 + 4Cz + (n, - 1) Cz] [1 - 4Cz - {ns - 1) Cz] 4 

= 4 . (3.26) 

Thus the variance has no dependence on the large-scale mode whatsoever. As a consequence, 
there can be no scale-dependent bias in single-field inflation^*^. 

We clearly see that this result is trivial. The transformation to fixed physical scale simply 
corresponds to evaluating the small-scale variance in the frame where there is no large-scale mode, 
i.e. it is the variance of (5s(x) on a scale R, which by definition does not depend on tpi. Thus, we 
could have written down the result that da'j^/dipi = 0 from the start. Instead, we had generated 
the mode coupling in a less convenient coordinate system and then transformed it back to the 
original coordinate system. 

Finally, we would like to again emphasize that these results are identical to those that we derive 
from perturbation theory in synchronous gauge, see Appendix B. Although the method used here 
(and in the next subsection) greatly simplify the calculation, all approaches give precisely the 
same result as long as one is careful to define the halos in terms of a fixed physical scale. 

3.3 Scale-dependent bias from direct application of consistency conditions 

A very useful observation emphasized in [31] (see also [30]) is that the procedure outlined above 
does not assume that the original background is homogeneous. Given a solution in a specific 
gauge, we can find a new solution to the equations of motion in the same gauge by performing 
a diffeomorphism. This procedure is linear in the long wavelength mode (and leading order 
in derivatives of the long mode), but is non-perturbative in the short wavelength modes. We 
used this feature in the previous subsections to determine the mixing between the long and 
short modes. However, there is no reason we cannot apply this procedure directly to the halo 
density field without ever discussing the biasing relative to the matter density. The halos may 
be determined by some non-perturbative behavior of the short modes, but to determine the 
influence of the long mode, all we need to know is how the halo density transforms under this 
diffeomorphism. 

One can simply repeat the argument from the pervious section replacing p by ph, the halo 
number density. Working in synchronous gauge, we again find the long mode via x = x (1 — (^;). 
The halo density transforms under a spatial diff as 

Ph{T) (1 -h 4(x, r)) = /9(x, r) = p(x, r) = p,,(r) (l + 4(x, r)^ , (3.27) 

so that 6s transforms as a scalar (as did Cs) under the dilation, 

4(x,r) = 4(x,r). (3.28) 

^°The Newtonian limit of this result was also explained in [36], in which case the coordinate transformation acts 
trivially. 
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Thus, in the coordinate system with the long mode explicit we find 

4(x, r) = 4(x(l + Cl)) = 4(x) + 0 X • V4(x) = 4(x) + 0{6hS). (3.29) 

Notice that the contribution from the long mode only enters at qnadratic order in the fluctuations. 
Therefore we can conclnde at linear order that 

(3.30) 

As discussed in Appendix A, in other gauges (e.g. Newtonian gauge) there may contributions 
at linear order that are proportional to drPh- However, these terms appear for densities at 
fixed coordinate time, rather than proper time. These linear contributions will cancel out when 
expressed in terms of observable quantities which are gauge invariant, as we explained in Section 2. 

4 Summary and Discussion 

Scale-dependent bias in the clustering of galaxies or halos is a powerful probe of primordial 
non-Gaussianity and therefore of the physics of inflation. Specifically, constraining local non- 
Gaussianity with a precision it(/nl) ^ 1 is, in principle, achievable with future surveys and may 
distinguish between single-field and multi-field inflation models. In order to correctly use such 
a measurement to test single-field inflation, an exact theoretical prescription for the observed 
galaxy overdensity, Ag, is needed. This prescription must take into account the fact that is 
measured in terms of observed redshifts and observed angular positions of galaxies, should be free 
of gauge ambiguities, and should include an exact prediction for the level of scale-dependent bias 
in single-field inflation. The latter question in particnlar has been a source of confusion and the 
main contribution of the present paper has been to clarify the calculation of this scale-dependent 
bias and to derive a consistent answer in several independent ways. The final resnlt, presented 
in Section 3, is that there is zero scale-dependent bias, Ab^a = 0, if inflation is governed by a 
single field, leading to the complete prescription for galaxy clustering snmmarized in Section 2. 

Most of the confusion abont scale-dependent bias in single-field models stems from the fact 
that matter pertnrbations in, for example, synchronons gauge have order-unity mode coupling, 
/nl = —5/3 — 5/12 [us — 1). Since scale-dependent bias arises from a modulation of the initial 
small-scale variance of matter perturbations (crj^) by a large-scale curvature perturbation (Ci), at 
first sight the mode coupling would thus appear to produce a scale-dependent bias corresponding 
to /nl of order unity. 

However, we have argued that this is not the correct interpretation. A crucial realization 
is that it can only be the initial variance of matter perturbations on a given physical scale, 
R = (1 + Cl) that determines the halo number density in that region at late times. While the 
initial variance of matter pertnrbations on a fixed coordinate scale does indeed see the order-unity 
modulation mentioned above, we have shown explicitly in Section 3.2 that when the variance is 
converted from being evalnated at fixed coordinate scale R to being evaluated at fixed physical 
scale R, the terms depending on Ci cancel and the modulation is identically zero (with even the 
term proportional to — 1 vanishing). 

To obtain the above result, we have derived the synchronous gauge mode coupling of matter 
perturbations in two ways. In Appendix B, we followed the approach common in the literature. 


esync. _ ?sync 
°h — °h,S 


A6ng(^) — 0 
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which matches the second order perturbation theory solution in the post-inflationary epoch to the 
inflationary consistency relation. However, in Section 3.2, we derived the same result in a simpler 
and, we argue, more insightful way. In single-field inflation, the local effect of a large-scale curva¬ 
ture perturbation Q (be. up to gradients) on the physics on small scales is equivalent to a simple 
coordinate transformation. In other words, there is no local physical effect. This follows from the 
fact that inflation is described by a single clock in the single-field case and each region of space 
finds itself on the same attractor solution. It remains true after inflation as long as perturbations 
evolve adiabatically, which is always true on super-horizon scales (more generally, it is true on 
scales above the sound horizon). Using this principle, which is equivalent to the inflationary con¬ 
sistency condition(s), we readily generated the well known mode coupling of synchronous gauge 
matter perturbations by a diffeomorphism, = —5/12din A|(A:)/dln = —5/3 — 5/12 (n* — 1) 
(see the derivation of Eq. (3.21)). In this picture, the aforementioned cancellation of the Q 
dependence of the variance became trivial: when converting to the initial variance of matter per¬ 
turbations on a fixed physical scale, we were simply undoing the diffeomorphism that produced 
the synchronous gauge mode coupling in the first place. 

Finally, in Section 3.3, we derived the absence of scale-dependent bias in an even simpler 
way by applying the diffeomorphism generating Q to the halo number density directly. We 
demonstrated that, modulo gradients, the long mode does not affect the halo number density at 
first order in perturbations and thus cannot generate a scale-dependent bias. 

We would like to emphasize that the absence of scale-dependent bias in single-field inflation 
has been to various degrees understood in the vast literature on consistency conditions in large- 
scale structure (e.g. [31-44]). It was for example appreciated by many authors that the effect of 
the long mode on various quantities is zero when expressed in terms of physical scales [57], in 
convenient coordinates like Fermi-normal coordinates [18, 44, 58] or in Lagrangian space [41]. The 
main contribution of the present paper was to resolve the paradox presented by the order unity 
mode coupling in coordinate space that led others to conclude there must be scale-dependent bias 
of order unity even in the single-field context. Moreover, we provided an explicit recipe for how 
to use this result to compute the observed galaxy perturbations taking into account relativistic 
terms (see also [42]). The simple recipe we provide can be applied in any gauge and can therefore 
be easily incorporated into any existing calculation. We thus now have the tools to properly 
interpret observational constraints on scale-dependent bias even when those constraints reach a 
precision (t(/nl) ^ 1- If such constraints favor non-zero scale-dependent bias as carefully defined 
in this paper, it would rule out the class of single-field inflation (modulo the caveats that we have 
assumed the inflationary evolution to have been on the attractor solution). 
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A Consistency Conditions in Newtonian and Synchronous Gauge 


In this appendix, we will review how the adiabatic modes arise from diffeomorphisms that leave 
Newtonian gauge fixed. We will then repeat the analysis in synchronous gauge and discuss the 
technical differences. Finally we will compare the description of halos in each gauge and their 
relation to physical observables. 

To derive the consistency conditions in conformal Newtonian gauge, we will follow [31] and 
define the metric to take the form 


ds‘^ = o^(f)[—(1 + + (1 “ 2^)6ijdx^dx^] . 

Now let us consider the transformation 

r = r + e(f) 

X = x{l — A) 


(A.l) 


(A.2) 

(A.3) 


where e(r) is an arbitrary function of time and A is a constant. This transformation shifts the 
potentials 


^{x, t) = l>(®, t) — e' — Tie 
'I’(®, r) = T(ic, r) — A + Ffe 


(A.4) 

(A.5) 


but leaves the gauge fixed. Under this change of coordinates, any scalar, s, will transform (to 
linear order) as 

s^(r, x) = s(t, x) — e(T) s'(t, x) , (A.6) 

where the is a reminder that this is the scalar fluctuation in Newtonian gauge and ' = dr- If 
we start from a solution to Einstein’s equations, any such transformation will generate a new 
solution. However, most such solutions cannot be extended to physical solutions with large, 
but finite wavelengths {e(r),A} —)■ {e(r, x), A(x)}. The reason is that some of the equations 
of motion may simply vanish when k = Q and are not satisfied for any A: 7 ^ 0. In conformal 
Newtonian gauge, the only such equation is — T) = 0. Therefore, if we impose = T as an 
additional constraint, this transformation can be extended to a physical solution. This requires 
that e' + 2Ffe = A or ^ 

e = ^ y dT'a^{T') = DyX , (A.7) 

where Dy is velocity growth function. In terms of C = —'k — H5p/p this means that C = C + -^ 
(where we used Sp = —p'e) and therefore we see that from a differmorphism we have generated a 
long wavelength Q. Therefore, if we start from a FRW solution and make the identification X = Q 
and C = 0 , we generate a new solution 


4 > = 
T = 
= 


Ci{D'^ + nDy) 

(A.8) 

Ciii-nDy) 

(A.9) 

Cl DyS' . 

(A.IO) 
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More generally, new solutions that are linear in Q can be generated from any solution for the 
short modes, even non-perturbative solutions. This procedure can also be extended to include a 
gradient of the long mode [31] (via the conformal consistency conditions [20, 55]) and to include 
contributions at non-linear order in the long mode [59, 60]. We can then use this procedure 
to determine the mode coupling between long and short modes provided only that we neglect 
contributions 0{didjCi). 

Although the details of the above argument were specific to Newtonian gauge, such an argu¬ 
ment will work in any gauge. For our purposes, we are are interested in synchronous gauge where 
the metric is taken to be 


ds^ = a^(r) —+ (1 “ 2^) 5^ dx^dx^ + didjxdx^dx^ 
It is easy to check that the transformation 


(A.ll) 


x = x{l-Ci) (A. 12 ) 

leaves the gauge unchanged and shifts the potential 

Ip ^ Ip- Cl ■ (A.13) 


Unlike Newtonian gauge, there is no change to the time coordinate and therefore the transfor¬ 
mation law for a scalar is simply 

ds® = 0 , (A.14) 

where ® is a reminder that this is the fluctuation in synchronous gauge. A priori, one may have 
been concerned that this was not a physical solution, as much of the gauge freedom has not been 
fixed by the choice of metric. Nevertheless, we can again check that in terms of C = —ip — H5p/p 
we again have C C + 0- This shows that the usual redundancies of synchronous gauge are 
playing no role in this transformation as those redundancies leave C fixed (after all, C is often 
called “gauge invariant” because it does not transform under the redundancies of synchronous 
gauge). One added benefit of this gauge redundancy in synchronous gauge is that we do not need 
to impose an extra condition on our transformation in order to extend it to a physical solution 
(one can check that there are no equations of motion that vanish when k = 0). 

A priori, one might think that to see theses results in different gauges, we may simply look 
up the diffeomorphism that takes us from one gauge to another. However, since our solution 
is equivalent to a diffeomorphism that leaves the gauge fixed, the diffeomorphism that takes us 
between gauges need not preserve the form of the long wavelength mode. As a result, to under¬ 
stand the implications of the long mode in a given gauge, it is most useful to work directly in the 
appropriate gauge from the beginning and add the long mode using the relevant diffeomorphism 
that keeps that gauge fixed. 

Now let us consider the implications for the matter density contrast <5 in synchronous gauge. 
The total matter density p{x) is a scalar under diffeomorphisms and therefore 

/(«, t) = p^six, t) + Ci{x)x • -k ... . (A.15) 
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We see that we have generated mode coupling between the short mode, p*, and the long wave¬ 
length mode, Cl, as described in the main text. We can perform a similar calculation in Newtonian 
gauge where we find 

P^{x,t) = p^{x,t) + CiDvP + Cl X- dp^ . (A. 16) 

Dehning ^ -1 where p is the homogenous solution, we find 

5^ = + Cl{x)x ■ d6l 5^ = 5^ + D,Ci{5' + ^{l + 5s)) + Cix-d6^ (A.17) 

In both gauges, we find a coupling between the long and short modes. However, in Newtonian 
gauge we hnd additional linear terms in Cl as well as extra mode coupling contributions which 
arise from the time component of the diffeormorphism. 

From the above procedure is should be clear that nothing required that p was the dark mat¬ 
ter density. It could have been any scalar quantity, including the density of halos or galaxies. 
Therefore, without doing any additional work, we see that at linear order in 5g and Cl we have 

= C + ^ A6 p,g = 0 (A.18) 

where d® is the physical density contrast. We emphasize that this means that, in the presence of 
a long wavelength mode, there is no change to density contrast of halos in synchronous gauge. 

One might worry that in Newtonian gauge, we do hnd terms linear in Cl that would appear 
as scale dependent bias. However, this feature arises because we are describing the density on 
constant coordinate time slices, rather than in terms of some physical dehnition of time. If we 
were to compute the observed galaxy density in a redshift survey, one hnds that the results appear 
in terms of the the gauge invariance density contrast 

(^invariant — (A.19) 

a In a " 

where the 5t arises from the fluctuations in the observed redshift relative to the coordinate time 
and therefore the last equality follows from the feature that dr® = 0. It is clear that there will 
be no scale dependent bias when written in terms of observable quantities in any gauge, but this 
is most transparent in synchronous gauge. 

B Second Order Perturbation Theory Calculation of Mode Coupling 

In Section 3.2, we used a diffeomorphism to directly generate the mode coupling of matter 
perturbations in synchronous gauge, = —5/3 — 5/12 {ug — 1). This result has been derived in 
the literature many times by solving the second order perturbation theory equations for the matter 
perturbations in the post-inflationary Universe and by then matching this to the inflationary 
consistency condition on C- Here, we briefly review this alternative derivation in the simple case 
of a matter-only Universe and show that it indeed leads to the same result, albeit through a more 
cumbersome calculation. We will not rederive the second order perturbation theory solution itself 
from scratch, but instead use well known expressions from the literature, specifically from [22, 61] 
(see also, e.g., [24, 62]). 
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Let us expand the perturbations to second order as <5 = (5^^^ + 1/2 5^^), X = X^^^ + l/2x^^\ 
'0 = ^(1) + 1/2 In synchronous gauge, in a matter-only universe, the first order solution is 
then given by 

t2 

5 

yVV (B.l) 

(note that we have fixed the remaining gauge freedom present in synchronous gauge). 

We split the solution into small- and large-scale components, (p = tps + ^pi- Then, to compute 
the mode coupling of interest, we need the second order part of the small-scale solution, but only 
those terms that are products of one small-scale mode and one large-scale mode, and only to 
zeroth order in gradients of (pi. Specifically, we care about terms of order tpi ~ ipir'^V'^cps 
contributing to and terms of order <piips contributing to The full calculation of the 

second order solution is cumbersome, if conceptually straightforward, and has been done for us 
(see, e.g., Eqs. (6) and (7) of [22]). Focusing on the relevant terms, the general result is the class 
of solutions^^ 

5(2) = y y 1 y v2^(2). (B.2) 

Note that, to the order of interest, simply solving the second order perturbation equations during 
matter domination allows us the freedom to choose 1 /( 2 ) (at least the relevant 0{(pi ips) contribu¬ 
tion), which is why 5 ( 2 ) is expressed in terms of it. For example, the equations would be solved 
by 

^ T’z y v2(^^ (B.3) 

^(2) = 0. 


^( 1 ) = 
^(b = 
5(b = 


However, we can always apply a redefinition, 

p)s = {1 +ipia)ips, (B.4) 

express Fq. (B.3) in terms of (ps and finally drop the tilde, to get a new second order solution, 

= Y y V2(avJs) 

= yT’zaT’s- 

Thus, indeed 5(2) = ‘fp>iT‘^ V^^Ps + V2^(2). 

While the matter domination second-order perturbation theory calculation thus does not 
uniquely specify the solution, we can fix it by requiring that the first order solution is Gaus¬ 
sian, which, in terms of mode coupling, means that is statistically independent of (note 

^^Note that [61] appears to have a typo causing an inconsistency between the solution for the metric perturbation 
and the matter overdensity. 
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that this is already done in the solution Eq. (6) and (7) in 
solution to the consistency relation for Q (see Section 3.1). 
perturbation on constant density slices is dehned as 

[22]). We do this by matching the 
Recall that the comoving curvature 


(B.5) 

We match the perturbation theory solution Eq. (B.2) to (/, in the regime where both the large 
and small modes are far outside the horizon. We can then ignore gradients acting on both ipi and 
^s- In this limit, (Sp=o = the right-hand side being the solution in synchronous gauge 

(this follows from, e.g., Eqs. (131) and (134) of [63]). Thus, 

F = -C-CF 

(B.6) 

or in terms of the mode coupling, 

^, = -C.(l + 2Cz). 

(B.7) 

The consistency relation for C, says 


C. = (i + Czx-V) cf. 

(B.8) 


Therefore, 

^, = -(l + 20 + 0x-V) Cf, (B.9) 

where CF is independent of Q so that this gives the true statistical mode coupling. Identifying 
the hrst order perturbation theory solution with the Gaussian held, i.e. = | (/j = —now 
hxes the second order solution to 




100 /, 1 
y (P, V^<P, - ^ 


20 2 ^2 10 2 


(l + ^ X • v) (/:1s 

+ ^X • V2(/9s 


• V(V2(/9s) 

_Hy,,^2x.v5W = 2Czx.V5F\ 


(B.IO) 


with (5^^ independent of the long mode by dehnition. 

We see that this is precisely the same result for the mode coupling that we derived directly 
using the consistency conditions, Eq. (3.19). As explained in Section 3.2, this is equivalent to 
local non-Gaussianity with = —5/3 — 5/12 {ug — 1). However, as detailed in Section 3.2, the 
mode coupling disappears when the short wavelength modes are evaluated on a hxed physical 
scale as opposed to a hxed (synchronous gauge) coordinate scale. 
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